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Non-linearity of the current-phase relationship of a Josephson junction is the key resource for a 
Josephson parametric amplifier (JPA), the only device in which the quantum limit has so far been 
achieved at microwave frequencies. A standard approach to describe JPA takes into account only 
the lowest order (cubic) non-linearity resulting in a DufRng-like oscillator equation of motion or in 
a Kerr-type non-linearity term in the Hamiltonian. In this paper we derive the quantum expression 
for the gain of JPA including all orders of the Josephson junction non-linearity in the linear response 
regime. We then analyse gain saturation effect for stronger signals within semi-classical approach. 
Our results reveal non-linear effects of higher orders and their implications for operation of a JPA. 


I. INTRODUCTION 

Non-linear nature of a Josephson junction has been 
recognized as a valuable resource for parametric amplifi¬ 
cation over several decades [1], The earliest experiments 
with a microwave oscillator incorporating a Josephson 
junction as a non-linear inductance have incontrovertibly 
shown that these devices, called Josephson parametric 
amplifiers can operate as amplifiers with extremely low 
noise temperatures [2] - [4]. Rapid development of low- 
temperature experiments with superconducting qubits 
requesting the quantum limited amplifiers to control and 
measure a broad range of parameters of quantum systems 
and absence of suitable commercial amplifiers stimulated 
great interest in JPA in the recent decade. Recent exper¬ 
iments confirmed that JPA can work as an amplifier near 
quantum limit [5] and can be used for a wide range of ap¬ 
plications. These include measurement of a displacement 
of a nano-mechanical oscillator [6], [7], implementation of 
the high-fidelity single-shot measurement [8]-[ll], quan¬ 
tum feedback [12], [13] and feedforward [14] for super¬ 
conducting qubits, as well as generation of the two-mode 
squeezed states of microwave photons [15]. 

A standard approach of a description of JPA takes 
into account only the leading order non-linearity origi¬ 
nated from a Josephson junction. In particular, keeping 
the first non-linear term only in the Josephson current- 
phase relationship one gets an oscillator equation of mo¬ 
tion with the cubic non-linearity analogous to the Duff¬ 
ing equation. In quantum approach the Hamiltonian 
acquires the Kerr-type non-linearity term coming from 
truncation of the energy of the Josephson junction. This 
simplification allows for an analytical solution for the 
steady state of an oscillator under strong pumping con¬ 
ditions. By linearising equation of motion for a weak 
external signal in the presence of a strong pump field one 
can obtain the analytical expressions for the gain of JPA 
(see for example [16]). These formulas provide reasonable 
agreement with experimental data but some discrepan¬ 
cies has been found, especially for low Q-factor oscillators 
[16]. In addition, the linear response theory did not ac¬ 
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count for gain saturation effects in the presence of strong 
signals [16], [17] which are essential for the understand¬ 
ing of limitations on the dynamic range of an amplifier, 
a property of large practical importance. 

The main goal of our paper is to reveal and analyse 
the non-linear effects in the JPA due to unavoidable pres¬ 
ence of the higher non-linear terms in the current-phase 
relationship of a Josephson junction. First, we account 
for all non-linear terms in the steady state solution of a 
driven oscillator at the pump frequency. This allows us 
to derive the exact quantum expressions for the gain of 
a parametric amplifier in the linear amplification regime. 
Secondly, we derive and solve the equation of motions be¬ 
yond the linear response theory to account for gain satu¬ 
ration effects under semi-classical approximation taking 
into account all orders of non-linearity. After revealing 
non-linear effects of higher orders we discuss applicabil¬ 
ity of the obtained solutions and their implications for 
operation of JPA. 

The paper is organized as follows. In Sec. II we write 
down the Hamiltonian of a JPA and formulate the input- 
output theory to derive the exact equation of motion. 
The response for a strong classical pump field account¬ 
ing for a complete sine non-linearity is found in Sec. III. 
In Sec. IV we derive equation for a weak quantum sig¬ 
nal and solve it analytically in the frequency domain. 
Having found the JPA gain expression we study higher 
order non-linear effects. In Sec. V we derive equation for 
an arbitrary strength quantum signal including both the 
fundamental non-linear terms and all non-linear terms 
and solve the equation numerically using semiclassical 
approach. The solution was used to describe the gain sat¬ 
uration effect in a JPA and estimate its dynamical range. 
Conclusions and final remarks are made in Sec. VI. 


II. EXACT HAMILTONIAN AND 
INPUT-OUTPUT RELATIONS 

The most basic JPA model consists of a parallel lossless 
LC resonant circuit where the Josephson junction plays 
a role of an inductor with a non-linear dependence on the 
current L(i) and the capacitance C comprises geometric 
capacitance of the circuit and the internal capacitance 
of the junction [Fig. 1(a)]. Without loss of generality 
we consider that the resonant circuit is galvanically con- 
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nected to a semi-infinite transmission line with character¬ 
istic impedance R (note that in practice the impedance of 
the transmission line is typically fixed to 50 Ohm and the 
effective impedance seen by the LC circuit is controlled 
by an additional coupling capacitor). 

With help of DC and AC Josephson effect equations, 
the Kirchhoff’s circuit laws and the Telegrapher’s equa¬ 
tions theory of a transmission line one can obtain an 
equation for the phase of the Josephson junction as 


d 2 ip 

dr 2 


1 dip 

o* +sm *’ 


i p cos (Dr) + i s (r), 


(1) 


where the amplitude of the harmonic pump current i p 
and the signal current i s (r) are normalized to the crit¬ 
ical current of the Josephson junction I c , Q = ujqRC 
is the quality factor of the oscillator, uiq = 1/y/LjC 
is the natural resonance frequency of the circuit, Lj = 
$ 0 /( 27 r/ c ) is the Josephson inductance and $ 0 = h/2e 
is the flux quantum. We also use dimensionless fre¬ 
quency Q = LOp/ujQ and time r = LOgt. By approximating 
sintp ~ ip — ip 3 /3 one obtains the equation of motion for 
Duffing oscillator used extensively to analyse parametric 
amplifier properties [17]. 
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FIG. 1: (Color online) (a). Circuit diagram of the most basic 
JPA. (b). Description of JPA with the input-output formal¬ 
ism. 


constant K defined as 


cuq 2 ewo 


( 5 ) 


Parametric amplification can be naturally studied in 
the framework of the input-output formalism [18], [19]. 
By introducing incoming Ai n and outgoing A out mi¬ 
crowave fields [Fig. 1(b)] the intra-oscillator A field is 
governed by the following quantum stochastic differen¬ 
tial equation 


dA 

dt 


A{t), H S ys 


\A(t ) + Vl A in (t), (6) 


where the coupling rate is given by 7 = wo/Q and we 
neglected the internal losses of the LC oscillator. The 
boundary conditions for all the fields at the input port 
of an oscillator read 


A out (t) = y/jA/t) - A in (t). (7) 

Substituting the Hamiltonian (4) into Eq (6) and tak¬ 
ing into account the expression for the commutator 

AAA + Ac = 2fc(A + ypJ for all integers 

k > 1 we obtain the following exact equation of motion 


dA 

dt 


= 0 X 


->2fc-l 


“[ (2 k + 1)! \wo 


K 


(i + A f ) 


2/c+l 


(^— + A(t) + 1/7 Ain(t). 


( 8 ) 


III. CLASSICAL RESPONSE AND 
STEADY-STATE SOLUTION FOR PUMP FIELD 


After a standard quantization procedure the quantum 
Hamiltonian of a capacitively shunted Josephson junction 
reads 

H = ^ +Ej{ l~cospl ( 2 ) 

where q is the operator of a charge on the capacitor and 
Ej = / c $ 0 /(27t) is the Josephson energy. 

By introducing creation and annihilation operators 


An input field of JPA consists of a strong harmonic 
pump field at the pump frequency uj p and a weak signal to 
be amplified. Extracting the harmonic time-dependence 
at the pump frequency uj p explicitly we respectively write 
the input, intra and output fields as 

A in (t) = (a in + a in (t)) e~ lUpt , 

A(t) = (a + a(t)) e ~ luJpt , (9) 

A ou t{t) = {aout + d ou t{t)) e~ lUpt . 



( 3 ) 


and dropping the constant zero-point motion term we 
convert (2) to 


H = 


tVjjQ 


A'A + Y. 


k—2 


22fc—3 

w 



( 4 ) 

Here we explicitly expressed the non-linearity as a sum of 
different orders of the weak effective Kerr non-linearity 


Here the first terms on the right hand sides of (9) repre¬ 
sent the strong pump fields treated as classical complex 
amplitudes (ai„, a, a ou t) while the other terms corre¬ 
spond to weak quantum signal fields described by oper¬ 
ators (d^, h, clout)' 

Neglecting weak signal contributions in (9) and substi¬ 
tuting only the classical pump fields into (8) we obtain 
the following equation 


— iuj p a = — lui^a 


E 

k=1 


~)2k—l 


k\(k + 1)! \ ujq 


K\c 


Q + iujo^j a + y/^a in 


( 10 ) 
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where we have taken into account the terms at the pump 
frequency only, applying the following approximation 

(ae~ iuipt + a*e iUpt ) 2k+1 + HU \a\ 2k e~ iuipt . 

v ' k\(k +1)! 1 1 

( 11 ) 

We note that the approximation (11) is equivalent to 
accounting only for terms with equal number of cre¬ 
ation and annihilation operators in the Hamiltonian (4). 
Eq. (10) has a harmonic steady-state solution at pump 
frequency uj p and is a good approximation of an exact 
solution for weak non-linear processes. 

It is well-known that the behaviour of the non-linear 
response (even for weak non-linearities) depends strongly 
on the amplitude of input field ai n [20], [21]. Specifi¬ 
cally, the resonance frequency deviates from the natural 
frequency ujq and decreases monotonically as the input 
field increases. At some critical value of the input pump 
field the response becomes bistable in some range of 
frequencies while below aZ lt the fields have unique and 
stable solutions at any frequency. The bistable regime 
(above critical pump) is used in Josephson bifurcation 
amplifiers [22] and is not considered in this manuscript. 

One may obtain the analytic expression for the crit¬ 
ical value of an input field ctf” 4 = ^/y/TJK (see 

for example [16]) if only the leading (quartic) non-linear 
term in the Hamiltonian (4) is kept. It is convenient to 
normalize the input field to this value as a.i n = retail™*, 
where r is the real scaling factor and <j> is the initial phase 
of the input pump field. Multiplying both sides of equa¬ 
tion (10) by their complex conjugate parts we obtain the 
equation 


1 

W + 


i-n + E 


fc=1 


2 2k ~ 1 (—n) k 
k\(k + 1)! 


2 ' 


n = 


VriQ 3 ' 


( 12 ) 

where H = uj p /loq is the dimensionless pump frequency 
and n = —|a| 2 AT/wo is the normalised number of pump 
photons in the oscillator. Eq. (12) has a convenient form 
if one needs to keep only a certain number of the non¬ 
linear orders in the series. For example, keeping only 
the terms for k = 1 we obtain the well-known result for 
the cubic non-linearity. One can use the following closed 
form of (12) to account for a complete sine non-linearity 


1 

4Q 2 



J\ (4y^n) \ 2 
4^ ) 


n = 


VriQ 3 ' 


(13) 


where J m (-) is the Bessel function of the first kind of 
order m. 

As mentioned above Eqs. (12)-(13) are normalized to 
the critical input field corresponding to the leading non¬ 
linear term. In other words, the case r = 1 corresponds 
to the bifurcation point for the cubic non-linearity only 
while the complete sine non-linearity has the bifurcation 
point at some other value of the parameter r. In order to 
find this value one can plot a stability diagram [Fig. 2] 
revealing parameter domains corresponding to a unique 
solution and to the bistable mode of Eq. (13). 



FIG. 2: (Color online) Frequency-amplitude (H, r) stability 
diagram based on numerical solution of Eq. (13) for a fixed 
value of the quality factor Q = 30. Yellow region corresponds 
to a unique stable solution while the green region shows pa¬ 
rameter range for three different solutions (two stable and one 
unstable solutions). The most right point of the green domain 
(0.96988, 1.0401) is the bifurcation point for the complete sine 
non-linearity. 


Having found the number of photons n we can express 
the amplitude of pump field a using equation (10) as 


1 

VQu o 1 
2 Q 

1 

y/Qio o 1 
2 Q 


re^aZ* 


1-0 


00 2 2fc_1 


(-»)* 


ti kKk + iy. 


re^aZ* 


O + 


Ji (4y / n) \ ' 

4v^ ) 


(14) 


In experiment the JPA can be characterized via mea¬ 
suring the reflection coefficient S\ \ = a ou t/otin■ Using 
the boundary conditions (7) and equation (14) we obtain 


Sn 


- iQ 


iQ 



J\ (4y / u) \ 
4y^ / 
J\ (4y / u) \ ' 

4 v^ ) 


(15) 


IV. THEORY OF LINEAR AMPLIFICATION 

As the powers of the input signal {a\ n ai n ) and am¬ 
plified signal ( a\ ut a ou t) are typically much smaller than 
the one of the input pump field lain] 2 , presence of a 
weak quantum signal does not affect classical solution 
(14). Known as the stiff pump approximation [23], this 
assumption allows us to derive the equation of motion for 
the signal in the linear response regime. 

One can substitute Eqs. (9) into Eq. (8) and keep only 
the terms linear in quantum signal a, <U. Applying both 
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approximation (11) for the pump signal and linearisation 
procedure for the quantum signal one can get 


/ ^ ^ , \ 2/c+l 

\A + AM = (ae~ luJpt + aV^*) 


n2fe+! 


Oh 

+ (2k + 1) (ae- 1 "'* + aV"**) (a(f)e 


+ a t (f)e iajpt ) + 0(a 2 ) —*• (2k + 1)! 

I „.| 2 fc „,2 


a a 


2/c 




r a(t) 


Q ,-| a |2(fe 1) 

(k- !)!(& + !)! 


fc!(jfc + l)! 

ot(t) ) e" iajpt . (16) 


Using Eq. (16) we obtain the following linearised equation 
of motion for the signal 


where the coefficients are defined as follows 


h 

h 



.a 2 K J 2 (4:y/ri) 

i --- 

c^o 2 n 


1 

2Q’ 


(17) 


(18) 

(19) 


Eq. (17) together with its Hermitian conjugate equa¬ 
tion forms a system of two linear differential equations 
with time-independent coefficients and can be solved in 
the frequency domain. Without loss of generality we 
consider the amplification process of a signal at some 
frequency co s and write the time dependence of intra res¬ 
onator field as superposition of their Fourier components 


a(r) 




a(A)e~ lAT dA, 


( 20 ) 


i r°° 

at(r) = -= / at(-A)e-^dA, (21) 

V 27T J — oo 


where A = (u> s — ui p )/coo is the normalized detuning of 
the signal from the pump frequency. Applying Fourier 
transform to Eq. (17) and its Hermitian conjugate pair 
we eventually arrive at the system of two linear algebraic 
equations with the solution 


i(A) = 


1 — (ll + iA)din(A) + hdl n (— A) 


y/Quj o (U + 4" — I2I2 


at(-A) = 


l I 2 djn (A) - (h + ?A)a- n (—A) 

y/QuJo (li 4- iA)(Z* -)- iA) — I 2 I 2 


( 22 ) 


Applying the boundary conditions (7) in the frequency 
domain we express the output signal Fourier component 
via the input one as 

a out (A) = g(A)din(A) + m(A)aJ n (-A), (23) 


where we defined 


3(A) = 


1 It+iA 

Q(h+iA)(lt + iA)-l 2 lZ 


m( A) 


1 I 2 

Q (li + iA)(l* + iA) — I 2 I 2 


(24) 


Coefficients g(A) and m(A) satisfy the relation G(A) = 
|<?(A)| 2 = |m(A)| 2 + 1 where G(A) can be identified as 
the gain of phase-insensitive linear amplifier [16]. 

In the traditional JPA theory the gain is calculated in 
cubic approximation by keeping only the first term k = 1 
in Eqs. (12) and (14). Now we have the exact expressions 
for the gain (24) where we account for all non-linear terms 
in the equation of motion (8). Using this result we can 
analyse the effect of higher order non-linearities on the 
gain of JPA. 







FIG. 3: (Color online) Gain as function of pump frequency for 
the fixed pump power r = 0.99, zero signal detuning A = 0 
and different quality factors Q of an oscillator. We also choose 
fo = 7 GHz and critical current 7 C = 2 gA. The colour 
represents different number of terms N taken into account in 
equation (8). N — 1 corresponds to the known case of qubic 
non-linearity, N = 2,3 account for the next two higher quintic 
and septic non-linearities, N = 00 represents the solution for 
a complete sin</> non-linearity. A/ indicates the frequently 
difference between the maximum of gain and AG shows the 
difference in maximal attainable gain for N = 1 and one for 
the exact solution N = 00 . 


In the Fig. 3 the gain of JPA at zero detuning G(0) 
is plotted as a function of the dimensionless pump fre¬ 
quency fl for the fixed pump power close to the critical 
value and different values of the Q-factor. For each value 
of Q-factor we take into account different number of non¬ 
linear terms N. To make a link to a real experimental 
situation we choose typical values for the natural fre¬ 
quency fo = wo/27t = 7 GHz and the critical current of 
the Josephson junction I c = 2 fiA and adjust the pump 
power to recover the experimentally relevant values of 
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maximal gains. 

Our results show that higher orders of non-linearity 
lead to a decrease and a frequency shift of the maximal 
achievable gain for a given pump strength and a qual¬ 
ity factor. As the quality factor decreases the steady 
state normalized photon number n corresponding to the 
critical pumping strength increases (see Eq. (13)) and 
one needs to take into account more terms in the expan¬ 
sion of the Bessel function (in physical terms it means 
that the current through the Josephson junction start 
to approach the critical current I c ). That explains our 
observation that the effect is more pronounced for lower 
values of Q and is practically negligible for Q > 150 (not 
shown). Within relevant parameter range the main con¬ 
tribution to the difference between the known solution 
for the fundamental cubic non-linearity (N = 1) and 
the exact solution is provided by the next higher quintic 
(TV = 2) non-linearity while the contributions from septic 
(TV = 3) and all other orders are not substantial. 





FIG. 4: (Color online) Signal gain dependences on pump fre¬ 
quency for the fixed Q = 10, zero detuning A = 0, natural 
frequency fo = 7 GHz and critical current I c = 2 /. iA. The 
colour represents different number of terms TV taken into ac¬ 
count in equation of (8). The pump amplitude r was ad¬ 
justed to achieve the same values of maximal gain. A f n with 
n = 1, 2, 3 indicate corresponding frequency shifts of the max¬ 
imal gain between solutions for TV = oo — 1; TV = oo — 2 and 
TV = oo — 3, respectively. 


For any value of Q-factor we can reduce the difference 
between maximal values of the signal gain by varying 
the level of applied pump power. In Fig. 4 one can see 
the gain dependences on pump frequency calculated for 
the fixed Q = 10, zero detuning A = 0 and different 
number of non-linear terms TV. For each of the solutions 
the pump amplitude was adjusted to achieve the same 
value of maximal gain while all other parameters were 
kept constant. Our results show that the shape of the 
curves are virtually identical and the only effect of the 
non-linear terms of higher orders is a shift of the reso¬ 
nance frequency (gain profile) of the JPA. 


In experiment it is hard to calibrate the exact value 
of power applied to a device in the absolute scale. It is 
also difficult to characterize the frequency shift due to the 
high-order non-linearities as the resonance frequency of 
an oscillator depends on the pumping power. This sug¬ 
gests that the contribution of higher-order non-linearities 
to gain of JPA in the linear amplification regime is not 
very significant and hard to observe experimentally. 

V. GAIN SATURATION EFFECT 

The linearised signal and idler gain expressions (24) 
are valid as long as the output power of amplified fields 
is much smaller than the pump field power. Firstly, we 
estimate gain saturation effect due to the fundamental 
cubic non-linearity. Using Eq. (8) with only one term in 
the sum for k = 1 and Eq. (9), where the classical pump 
field a is already found via Eqs. (13)-(14), we obtain the 
following equation of motion for quantum signal in the 
cubic approximation 


— = l id + ha* + ci a 2 + C 2 (da^ + d^a) 

CLT 


+ c 3 (a 2 a^ + aa*a + a*a 2 ) + ^ . 

vQw o 

(25) 

The linear coefficients l\ and I 2 are given in 
while the cubic coefficients are defined as 

(18) and (19) 

K Ji(Ay/n) 

(26) 

ci - 1 a , 

w 0 2 y/n 

a 

c 2 = —Cl, 

cr 

(27) 

K J 0 (4y/n) 

C 3 = l q 

(28) 

cjo 3 


We note that to derive Eq. (25), an approximation sim¬ 
ilar to (16) is used again, where both the linear terms 
of quantum signal and the three non-linear terms cor¬ 
responding to the cubic non-linearity are taken into ac¬ 
count. In this way one can derive more precise equations 
for the quantum signal containing the non-linear terms 
of higher orders. Also, we can include all the non-linear 
terms in Eq. (17) to relax the restriction on smallness am¬ 
plitude of input signal. However, in this case the signal 
equation becomes cumbersome and not portable. For¬ 
tunately, we can simplify the signal equation conducting 
the derivation of the signal equation in the semi-classical 
limit, i.e. we take a(t ) to be a complex function and 
not an operator. Mathematically, it results in the two 
following approximations 

(A + Atj —► [(a + a(i)) e" 1 "’* 

+ (a* +a*{t))e iul ^] 2k+1 

-A (a + a(t)) k+1 (a* + e (29) 

On the first step we go to the semi-classical limit while on 
the second step we account for the terms ar the pump fre¬ 
quency only. Applying the approximations (29) to each 
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of the non-linear terms in Eq. (8) we finally get the exact 
non-linear signal equation in the semi-classical limit as 


da 

(It 


1 

2 Q 


1 r> 

2 4 sjm{a) J 


a{r) 


Ji(4y/n) Ji(4^/ro(a)) a in {r) 

—:— i= - , — a-f- . = (oU) 

4 yfn 4y/ m(a) y/Quio 


where m(a) = n — K (a*a + aa* + a* a )/wq ■ 

We used the fourth-order finite-difference (Runge- 
Kutta) method to solve the system of Eq. (30) and its 
complex conjugate pair as well as the system of the ap¬ 
proximate equation of motion (25) and its Hermitian con¬ 
jugate pair in the semi-classical limit. To calculate the 
gain of JPA we choose the excitation on the right hand 
side of Eq. (30) [or Eq. (25)] in the form of a harmonic 
signal at some frequency ui s : ai n {r) = ai n (A)e" iAr while 
in the complex conjugate pair of Eq. (30) [or (25)] we set 
a\ n (— A) = 0. By numerically calculating a ou t at the 
same frequency w s we obtain a ou t{ A) and evaluate the 
gain using Eq. (23). 

In Fig. 5 we present the maximal signal gain maxG(O) 
as function of the input signal amplitude ai n ( 0) corre¬ 
sponding to the zero frequency detuning A = 0 between 
the signal w s and the pump oj p frequencies. We use sim¬ 
ilar parameters of JPA as in the middle panel of Fig. 3. 
The red curve represents the solution of exact equation 
of motion (30) while the green curve correspond to the 
solution of Eq. (25) in the semi-classical approximation. 
We also solve Eq. (25) with C3 term only (blue curve in 
Fig. 5) by setting manually C\= ci = 0. For the solution 
of (25) we set the pump power close to the bifurcation 
point putting r = 0.99. In order to achieve the same 
gain level in the case when all non-linear terms (N = 00 ) 
are taken into consideration in the signal equation (30) 
we slightly change the pump power to r = 1.00297 (see 
discussion in the previous section). The analytical solu¬ 
tion (24) corresponding to the gain in linear amplification 
regime is shown as a horizontal line for reference. 

As expected, Fig. 5 demonstrates that the gain de¬ 
creases at large powers of input signals. The verti¬ 
cal dotted line in Fig. 5 corresponds to the relation 
ppump/psignni _ 2 q which represents the practically 

found rule for validity of the stiff pump approximation. 
We can also conclude that the fundamental contribu¬ 
tion comes from the third term in the cubic non-linearity 
3 c3 a*a 2 (blue line in Fig. 5). 

It is logical to expect that the dynamic range increases 
with decreasing the coefficients ci, C2, C3. Using the for¬ 
mula (28), expressing the critical current as I c = uiQCjipo, 
taking into account (5) and the quality factor formula 
Q = uJoCjR we find Cj oc K/ujq ~ l/Q. In order to il¬ 
lustrate this dependence we plot the maximal gain for 
three different values of the ratio ujq/KQ = —1, —10, 
— 100 [Fig. 6] by adjusting the critical current values keep¬ 
ing other parameters constant. It is evident that the 
JPA dynamical range increases with increasing the ratio 

u>o/\K\Q. 

Our solution is based on approximation that both the 



fflxn(O), s ! 


FIG. 5: (Color online) Signal gain saturation due to sine and 
cubic non-linearities. Here the natural frequency fo=7 GHz, 
the critical current I c = 2 ptA and the quality factor Q = 30 
of the JPA are the same as in the linear case depicted in the 
middle panel of Fig. 3. 



FIG. 6: (Color online) Dependence of the dynamic range on 
the Kerr non-linearity under constant the natural frequency 
fo = 7 GHz and the quality factor Q = 30. 


input pump field and the intra resonator (output) pump 
field have the same harmonic dependence in time (see 
Eq. (11)) which is justified for i p = I p /I c « 1. In¬ 
deed, equation (1) always has a period-one (harmonic) 
steady state solution at the pump frequency if the con¬ 
dition of small oscillations I p /I c « 1 is satisfied [21], 
[24]. With increasing the ratio I p /I c the form of solu¬ 
tion to Eq. ((1)) can be changed dramatically demon¬ 
strating a period doubling cascade through a sequence of 
pitchfork bifurcations (transition from regular to chaotic 
dynamics) and eventually chaotic behaviour [21]. Under 
these conditions the classical solution for intra resonator 
field (14) fails independently of the number of non-linear 
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terms taken into account. The number of pump photons 
in the JPA resonator at the point close to bifurcation 
is proportional to N p ~ a^ rit ~ IJQ 1 2 3 4 5 6 7 8 9 10 11 12 - Therefore, the 
pump current approaches the critical current when suffi¬ 
ciently low quality factor is chosen affecting the expected 
operation of the amplifier and imposing limitation on its 
dynamic range and bandwidth [16]. 

VI. CONCLUSION 

We obtained the expression for the gain of JPA includ¬ 
ing all orders of non-linearity in the linear amplification 
regime. Comparing to the conventional solution for the 
qubic non-linearity we found that higher order terms lead 
to a decrease of the maximum achievable gain for given 
pump power and to a frequency shift of maximum gain 
point. By increasing the pump power one can restore the 
value of the maximal gain and only the frequency shift 
of the maximal gain point remains. 

To study the effect of the gain saturation at high sig¬ 
nal powers we solved the equation of motion including 
all non-linearities in the semi-classical approximation be¬ 
yond linear response regime. We revealed the principal 
term responsible for the gain saturation effect and identi¬ 
fied the parameter controlling the dynamic range of JPA. 

Our solution is valid when both the input pump fields 
inside and outside of the resonator have the same har¬ 
monic dependence which is the case when the current 


though the Josephson junction is much less than the crit¬ 
ical current I p /I c « 1. In this case of weak non-linear 
process the equation of motion always has a period-one 
(harmonic) steady-state solution at the pump frequency 
and our approximation is justified. 

At higher driving conditions the solution to Eq (1) 
may demonstrate a period doubling cascade through a 
sequence of pitchfork bifurcations (like it takes place for 
the Duffing equation [21]) and our approach (as well as 
the conventional cubic formula) will be inaccurate. It is 
interesting to note that at higher driving conditions the 
solution to the Duffing equation has a parameter domain 
when a large amplitude period-one steady state solution 
[21]. Therefore, Eq. (1) can in principle also have the 
same properties. We speculate that in this case condi¬ 
tions for JPA operation can be found even for I p < I c 
which will allow simultaneous large bandwidth and large 
dynamic range. This will be a subject of future investi¬ 
gation. 
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